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Abstract. - The high-temperature aspects of the Casimir force between two neutral conduct- 
ing walls are studied. The mathematical model of "inert" ideal-conductor walls, considered 
in the original formulations of the Casimir effect, is based on the universal properties of the 
electromagnetic radiation in the vacuum between the conductors, with zero boundary condi- 
tions for the tangential components of the electric field on the walls. This formulation seems 
to be in agreement with experiments on metallic conductors at room temperature. At high 
temperatures or large distances, at least, fluctuations of the electric field are present in the 
bulk and at the surface of a particle system forming the walls, even in the high-density limit: 
"living" ideal conductors. This makes the enforcement of the inert boundary conditions inade- 
quate. Within a hierarchy of length scales, the high-temperature Casimir force is shown to be 
entirely determined by the thermal fluctuations in the conducting walls, modelled microscopi- 
cally by classical Coulomb fluids in the Debye-Hiickel regime. The semi-classical regime, in the 
framework of quantum electrodynamics is studied in the companion letter by P.R.Buenzli and 
Ph. A. Martin [1]. 



This letter is related to the one by Buenzli and Martin [1]. For the sake of completeness, 
we cannot avoid repeating a few things. 

Casimir showed in his famous paper [2] that fluctuations of the electromagnetic field in 
vacuum can be detected and quantitatively estimated via the measurement of a macroscopic 
attractive force between two parallel neutral metallic plates; for a nice introduction to the 
Casimir effect see [3] and for an exhaustive review see [4]. 

Let us recall briefly, within the formalism of Ref. [3] , some aspects of the usual theory for 
plates considered as made of ideal conductors, which are relevant in view of the present letter. 
We consider the 3D Cartesian space of points r = [x, y, z) where a vacuum is localized in the 
subspace A = {r\x 6 {—d/2, d/2); (y, z) £ R 2 } between two ideal-conductor walls (thick slabs) 
at a distance d from each other. The time-dependent electric E(r,i) and magnetic B(r,t) 
fields in A are the solutions of the Maxwell equations in vacuum, subject to the boundary 
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conditions that the tangential components of the electric field vanish at the ideal-conductor 
walls OA = {r\x = ±d/2; (y, z) e R 2 }: 

E y {r,t) = E z {r,t) = for r e dA. (1) 

Note that this mathematical definition of the ideal-conductor wall is based on macroscopic 
electrostatics: the electric field is considered to be zero, without any fluctuation, inside the 
walls which have no microscopic structure and act only as fixing the instantaneous boundary 
conditions of type QJ. We shall call such a mathematical model of an ideal conductor "inert 
ideal-conductor". For each separate mode labeled by the wave number k = (k x , k y ,k z ) with 
k x = im x /d (n x = 0,1,2,...) and polarization indices A = 1,2 (only one polarization is 
possible when k x =0), the quantized energy spectrum of the electromagnetic field between 
the walls corresponds to that of an oscillator with the frequency Uf, — c\k\ (c is the velocity 
of light). At zero temperature T = 0, no photons are present and so each mode contributes 
by the zero-point energy fuvk/2 where h is Planck's constant. The d-dependent part of the 
system ground-state energy leads to the following attractive Casimir force per unit surface of 
one of the walls 

At nonzero temperature T > 0, all numbers of photons are possible and each mode con- 
tributes by the free energy of the thermalized harmonic oscillator. The Casimir force then 
reads 

9 °° r°° 

/ T ( d ) = -iE'/ M ± k ± q n (e 2 ^-iy X , (3) 

" n=0 







where f3 = l/(/cgT) is the inverse temperature, the prime in the sum over n = 0, 1,2, . . . 
means that the n = term should be multiplied by 1/2, k± is the magnitude of a wave- vector 
component in the (y, z)-plane and — k\ +in/° 2 with £„ = 2im/(hf3) being the Matsubara 
frequencies. By a simple change of variables, formula © can be rewritten as follows 

where 

Aird , . 

is the dimensionless parameter which measures the ratio of the separation between the con- 
ductor walls to the thermal wavelength of a photon. The small values of t correspond to low 
temperatures or small distances where quantum effects dominate. Using the Euler-MacLaurin 
sum formula, one obtains from eq. (@J) the small-i expansion of the form 

It is interesting that the leading correction to the T = result ((2J) is negligible in the ex- 
periments which have been performed at room temperature, see for example refs. [5,6]. The 
experiments are in good agreement with JJjJ . The large values of t correspond to high temper- 
atures or large distances where the classical limit of quantum mechanics provides an adequate 
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system description. In the large-t limit, the n = term dominates in the sum which 
implies the classical ^-independent leading behavior 



C(3) 



1 



An(3d 3 13d- 



t 



(7) 



For the present time, the high-t region is not accessible to experiments on metals. However, 
the high-temperature regime might be of interest for electrolytes. 

Lifshitz [7] considered the more general case of dielectric walls with a frequency-dependent 
dielectric permittivity e(uj). His starting point was the fluctuations within the walls, which 
therefore were not considered as inert. He derived the following formula for the Casimir 
force [4]: 



H n=0 



dk± k±q rl 



- 1 



-1 



(8) 



where and r± are the reflection coefficients of the TM and TE modes, respectively. They 
are given by 



»"|| 2 (£n,k±) 



e(i£n)<?n 



e(i£n)<7« 



r ± 2 (£n,kj_) 



In 



(9) 



with fc 2 = k\ + e(i£„)£ 2 /c 2 . When e(uj) < oo, eqs. (JSJ and © are well defined. When 
e(ui) — ► oo, the zero-frequency n = term in the sum on the rhs of JSJ) is not uniquely defined 
because its value depends on the order of the limits e(u>) — ► oo and n — > 0. In order to 
restore the inert ideal-conductor result based on the electrostatic boundary conditions 
CJ, Schwinger et al. [8] postulated the following order: set first e(u) — oo, then take the 
limit n = 0. This prescription implies the reflection coefficients of the zero mode to be 
(0, k±) = r^(0, k±) = 1 for inert ideal metals. 
Experiments are performed on real conductors composed of quantum particles, with finite 
static conductivity a and plasma frequency oj p , given by oj^ = 47re 2 n/m where n is the number 



density of free electrons of mass m. 
the frequency-dependent e(w): 



For such real conductors, one has the Drude formulae for 
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for uo — > , 



e(w) - 1 1 for cj > cj 2 /(47rcr) . 

or ^ 



(10) 
(11) 



The consideration of a frequency-dependent e(a;) enables one to avoid an artificial prescription 
for the order of limits: it is the dynamics of the particle system which "chooses" the correct 
treatment of the zero-mode contribution. In a series of recent works [9-13], the Drude formula 
lfTT)|l was substituted into eq. © considered for the zero Matsubara frequency £o — * 0. This 
leads to the reflection coefficients f 2 (0, k±) — 1, r^(0, k±) = independent of a, i.e. for n = 
the second term on the rhs of eq. © does not contribute to the Casimir force for a real 
conductor. As a mathematical consequence, the additional term ((3)/(8irf3d 3 ) appears in the 
Casimir force in any regime. In particular, the large-temperature formula is modified to 



(12) 
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a result identical to the one given by Lifshitz [7]. Although the additional term vanishes at 
zero temperature, it is relevant in the region of small temperatures where it is the source 
of some contradictions. Namely, it was argued in another series of works [14-16] that, at 
low temperatures, the relation 1)11(1 should be used. An intensive polemic about the low- 
temperature Casimir effect persists in our days [13,16]. 

In this letter, we shall concentrate on the high-temperature aspects of the Casimir effect. 
There is an apparent discrepancy by a factor 1/2 between the high-temperature Schwinger 
formula (J7J, valid for inert ideal-conductor walls with the boundary conditions QJ, and the 
Lifshitz formula 112(1 . valid for real-conductor walls with e(w) given by the Drude dispersion 
relation (|1L)|) . We aim at explaining this discrepancy on the basis of some exact results 
for specific microscopic particle systems which are used to model the conductor walls. The 
consideration of the Casimir effect in the t — > oo limit is also motivated by two fundamental 
simplifications of these model systems. First, according to the correspondence principle, in a 
microscopic model of matter coupled to electromagnetic radiation at equilibrium, both matter 
and radiation can be treated classically in the high-temperature limit. This fact manifests 
itself as the absence of Ti in the leading terms of the expansions (0 and Q12(l . Second, the 
application of the Bohr- van Leeuwen theorem [17,18] leads to the decoupling between classical 
matter and radiation, and to an effective elimination of the magnetic forces in the matter (for 
a nice detailed treatment of this subject, see ref. [19]). The absence of relativistic effects is 
seen via the independence of the leading terms in eqs. (0) and ((12(1 on c. We conclude that 
the matter can be treated in the f -> oo limit as a classical matter, unaffected by radiation, 
where the charges interact only via the instantaneous Coulomb potential. 

As a model system of the classical Coulomb fluid, we consider a general mixture of M 
species of mobile pointlike (structureless) particles a = 1, 2, . . . with the corresponding masses 
m a and charges Z a e, where e is the elementary charge and Z denotes integer valence (Z = 
— 1 for an electron). Its statistical mechanics is treated in the grand canonical ensemble 
characterized by the inverse temperature and by the species fugacities {z a } or, equivalently, 
the bulk species densities {n a } constrained by the neutrality condition ^ Q Z a n a — 0. The 
thermal average will be denoted by (•••). We use Gaussian units. The interaction energy of 
particles {i} with charges {<?i}, localized at spatial positions {r^}, is given by J2i<j[ c li < lj v (\ r i ~~ 
rj\) + U(A i +A J )/2(l r i ~~ r-j ] )] , where v(r) — 1/r is the Coulomb potential and 



is the hard-core repulsion potential which prevents the classical thermodynamic collapse be- 
tween oppositely charged particles. To make the correspondence with the quantum-mechanical 
version of the model, the hard-core diameter of particles of type a has to be set equal to the 
thermal de Broglie wavelength A Q = h(27r(3/m a ) 1 / 2 [20]. 

We would like to emphasize that the present particle system represents a microscopic model 
of "living conductors" where the charge density and the corresponding electric potential/field 
fluctuate, even for extreme values of physical parameters like the particle density. To be more 
precise, let us consider the truncated charge-charge correlation function 



where the microscopic charge density p is defined by p{r) = ^ a Z a eh a (r) with n a (r) = 
S(Z a e, qi)5(r — r,) being the microscopic number density of a-species. In the infinite 
space (bulk regime) , the fact that the Fourier transform of the Coulomb interaction has the 




for r < A , 
for r > A , 



(13) 



S(r,r') = (p(r)p(r')) T 



(14) 
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form v(k) = Air/k 2 implies the following small-fc behavior of the charge structure function 
(the Fourier transform of (|14|) with respect to \r — r'|) 

S( k ) = ^ k2 + 0(k 4 ); (15) 

for a review of sum rules for charged systems, see ref. [21]. This exact result is not influenced 
by short-range interaction potentials like the hard-core one. Thus, the second moment of S(r) 
does not depend on the total particle number density n — ^2 a {fi a (r)) and survives also in the 
high-density region. An immediate consequence of eq. H15(l is an asymptotic formula for the 
long-ranged potential-potential correlation function [22] : 

/3(</>(r)0(r')) T ~ -, 7T as|r-r'|^oo. (16) 

\r — r'\ 

Here, (j>(r) — J dr'u(|r — r'\)p(r') is the microscopic electric potential created at point r 
by the system of charged particles, and the distance \r — r'\ has to be large compared to 
the microscopic scale represented by the correlation length of the short-ranged (exponentially 
decaying) particle correlations. Since the microscopic electric field E is given by E^(r) = 
—d^ir) (p = x,y,z), the field-field correlation function is obtained from eq. (|16|l as 

P(Mr)Mr')f ~ 3(r ' r%(r |r r,) ; | ;^' r ' r12 . (17) 

It is obvious that (E(r)} = since the mean electric potential is a constant inside a conductor. 
However, the asymptotic formula (|17|l tells us that nonzero thermal fluctuations of the electric 
field must be present in the system for any particle density n. The generalization of the 
fluctuation results, obtained for the bulk, to inhomogeneous situations of the present interest, 
like conductors with boundaries, was made in ref. [23]. As soon as the two points r and r' 
are inside a conductor, asymptotic formulae Ijl6|l and l|17(l remain valid. When one of the 
points lies on the conductor boundary, the tangential components of the electric field at this 
point still fluctuate according to l|17|) , while the discontinuity of the normal component across 
the surface is related to surface charge fluctuations. These fluctuation phenomena make the 
living conductors fundamentally different from the inert ones with tangential components of 
the electric field at a boundary identically set to zero, as in eq. Q). The Casimir force <|12|) 
can be retrieved through a Maxwell stress tensor computed from the electric-field fluctuations 
in the vacuum region [24]. 

As was already mentioned, our Coulomb fluid of classical charged particles with de Broglie 
hard cores can represent its quantum counterpart of pointlike charges at sufficiently large 
temperatures. It has been shown in [1] that the long-range charge correlations of the semi- 
classical regime do not spoil the classical limit l|12l) . The high-temperature region of classical 
fluids is described exactly by the Debye-Huckel (DH) theory. Rigorous conditions, under which 
the DH approximation gives the exact leading correction to the ideal gas, were the subject of 
many studies in the past; for a short historical review, see e.g. [20]. These conditions arise 
naturally in a renormalized Mayer diagrammatic expansion for statistical quantities [25,26]. 
In terms of the mean interparticle distance a and the inverse Debye length k (k _1 is the 
correlation length of particles in the DH regime), defined by 

^ = -, « 2 =We 2 ^Z> a , (18) 
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the DH scaling regime is given by (see Eq. (11) of ref. [26]) 

Q) «^ e2<<1 ' 

where A represents the "typical", in our case dc Broglie, hard-core radius of particles. Suppos- 
ing in what follows for the sake of simplicity that ^ Q Z 2 n a /n is of order of unity and omitting 
irrelevant numerical factors, these inequalities can be rewritten in a more transparent form 

a <C a < kT 1 (20) 

where eto ~ Ti 2 /(me 2 ) with m being the "typical" particle mass. The lightest of the charges 
are the electrons for which the quantum microscopic scale do attains its maximum value, equal 
to the Bohr radius ~ 10~ 10 m. The first "Bohr" inequality in l|2U|l is a quantum upper bound 
for possible values of particle densities. Since the Bohr radius is small, very dense Coulomb 
fluids with a ~ 10~ 8 m are allowed; we shall refer to them as "living ideal conductors". The 
second inequality combines both the particle density n and the temperature parameter /3e 2 ; 
for a fixed particle density allowed by the first inequality, there always exists a sufficiently 
high temperature above which this inequality is fulfilled. 

Let our classical Coulomb fluid model the conductor slabs in the Casimir geometry. The 
characteristic correlation length of the particle system is assumed to be much smaller than 
the macroscopic vacuum distance d between the conductor walls, 

K~ x <d. (21) 

Then, the i-parameter (0 is of the form A(nd)/(Ka) 3 with A of the order of 12Tre 2 / (Tic) , i.e. 
unity. The scale hierarchy l|20ll and (|21|l is thus fully consistent with the classical and nonrel- 
ativistic limit of interest t — > oo. To summarize: as soon as the scaling length regimes <|2(J[) 
and (|21|l apply, the quantum system of charged particles coupled to electromagnetic radiation 
at equilibrium can be represented in terms of its classical pure-Coulomb fluid counterpart, 
decoupled from radiation and treated within the DH theory. The Casimir force originates 
exclusively from the thermal fluctuations in the conducting walls modelled by this classical 
Coulomb fluid. 

The Casimir problem of microscopic Coulomb fluids was solved by using an inhomogeneous 
version of the DH theory in two recent papers: the work [24] dealt also with more complex 
physical situations, the study [27] went beyond the DH theory. In the DH theory, the large-red 
expansion of the Casimir force was obtained in the form [24] 

The leading universal term is identical to the Lifshitz result 1)12(1. The subleading correction 
term is non-universal and depends on the composition of the Coulomb fluid via k. Even for 
a very dense Coulomb fluid with the mean interparticle distance a ~ 10 _8 m (living ideal 
conductor), there exists a sufficiently high temperature and a sufficiently large distance be- 
tween slabs above which the required length scale hierarchy a <C kT x <rf takes place and the 
correction term is negligibly small in comparison with the leading one. 

In conclusion, the mathematical model of inert ideal-conductor walls is based on the zero 
boundary conditions for the tangential components of the electric field QJ. This seems to be 
in agreement with experimental results at zero temperature, and perhaps also at sufficiently 
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small temperatures. Why the quantum ground-state fluctuations in the real walls seem to play 
no role is an open problem. At high temperatures, fluctuations of the electric field prevail 
in the bulk and at the surface of the particle system, even in the high-density limit (living 
ideal conductor), which makes the enforcement of the inert boundary conditions inadequate. 
Within the hierarchy of length scales 1)20(1 and the high-temperature Casimir force was 
shown to be entirely determined by the thermal fluctuations of the conducting walls, modelled 
microscopically by classical Coulomb fluids in the Debye-Hiickel regime. 

* * * 
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